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Abstract 

The transition in quantum chromodynamics (QCD) from hadronic matter to the quark-gluon plasma 
(QGP) at high temperatures and/or net-baryon densities is associated with the restoration of chiral 
symmetry and can be investigated in the laboratory via heavy-ion collisions. We study this chiral 
transition within the functional renormalization group (FRG) approach applied to the two-flavor 
version of the extended Linear Sigma Model (eLSM). The eLSM is an effective model for the strong 
interaction and features besides scalar and pseudoscalar degrees of freedom also vector and axial- 
vector mesons. We discuss the impact of the quark masses and the axial anomaly on the order of the 
chiral transition. We also confirm the degeneracy of the masses of chiral partners above the transition 
temperature. We find that the mass of the ar meson (p meson) decreases (increases) towards the 
chiral transition. 


1 Introduction 

QGD is the fundamental theory of the strong interaction. For massless quarks, the QCD Lagrangian 
has a global U{Nf)ji x U^Nf)^ = SU{Nf)v x SU{Nf)A x U{l)v x U{1)a symmetry, where Nf denotes 
the number of quark flavors. The U{l)v symmetry corresponds to baryon-number conservation, which is 
always respected. At the quantum level, U{1)a is broken to Z(iVy)^, a phenomenon which is referred to 
as the axial anomaly [ij. In the QCD vacuum, the remaining SU{Nf)v x SU{Nf)A symmetry, termed 
“chiral symmetry” in the following, is further spontaneously broken to SU{Nf)v by a non-vanishing 
quark condensate (qq), inducing Nj — 1 Goldstone bosons |2|, [sl Ij] . For nonzero and degenerate quark 
masses, the chiral symmetry is also explicitly broken to SU{Nf)v- 

At high temperatures and/or net-baryon number densities, the quark condensate melts and chiral 
symmetry is effectively restored. This chiral transition is commonly associated with the so-called QCD 
transition between a hadronic phase and the QGP. The QGP state has existed during the early stages 
of the universe. Experiments at accelerator facilities, such as SPS and LHC at CERN, RHIC at BNL, 
or SIS-100/300 at the FAIR project in Darmstadt, aim to explore the QGP via heavy-ion collisions Q. 
Above the chiral transition, the masses of chiral partners, such as the sigma and the pion or the p and 
the oi meson, become degenerate Q. In particular, dropping p and ai meson masses were suggested as 
signatures for chiral symmetry restoration 0,0,0 • The change of the spectral properties of the p meson 
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during the chiral transition could, e.g., be detected via its decay into dileptons [1^. Modifications of the 


dilepton spectrum have been observed in Pb+Pb [ll| and In+In [12| collisions. Concerning the axial 
anomaly, it was claimed in Ref. [l^ that data for Au+Au collisions at RHIC energies show a reduction 
of the r]' meson mass, which was interpreted as a precursor for an effective restoration of the 11(1)a 
symmetry [l3| . 

In general, the order of the chiral transition depends on the number of quark flavors and their 
masses Il5l| . Furthermore, for Nf = 2, the order depends on whether the U{1)a symmetry is restored 


prior to the restoration of chiral symmetry llal. o r whether it remains broken by the axial anomaly. For 
vanishing quark masses, Pisarski and Wilczek [l7| have argued that, in the absence of the U{1)a anomaly, 
the chiral phase transition for Nf > 2 is of first order and, in the presence of the U{1)a anomaly, can 
be of second order for Nf = 2. For an infinite anomaly strength and if the transition is of second order, 
it definitely falls into the 0(4)-universality class [as originally conjectured by the authors of Ref. A 
Recently, however, it was argued that the chiral transition can be of second order for Nf = 2, no matter 
whether the axial anomaly is present 181 or not [l^, [l^, only the associated universality class differs: it 
was conjectured that the second-order transition can be in the SU{2)a x U{2)v class for finite anomaly 
strength, and in the U{2)a x U{2)v class for zero anomaly strength, respectively. For nonzero quark 
masses, the second-order transition is smeared out into a crossover. 

The chiral transition has been extensively investigated from first principles via lattice QCD 22, 
24, 25, 26, 27, 2^ 2^. The critical temperature for vanishing quark-chemical potential, fJ. = 0, was 
160 MeV 0, [s^, im, [H, 1^ . Lattice-QCD calculations do n ot agree upon whether the 

[H, 113, [ail. 


estimated to be 

U{\)a symmetry has already been restored at the chiral transition temperature 


or not 


Unfortunately, perturbative methods fail to achieve reliable results in the low-energy range, where the 
gauge coupling is strong. Besides, they are plagued by severe infrared divergences. Therefore, non- 
perturbative continuum methods, such as the functional renormalization group (FRG) [H, Ei, El, Ell, 
are complementary to lattice QCD and hence indispensable for exploring the nature of QCD matter. 

FRC studies of the chiral transition using effective models for QCD have a long tradition [4|: a 
quark-meson model with U(2)r x U(2)l symmetry has been investigated beyond the local potential 
approximation (LPA) in Ref. [dj,E|- An 0(iV)-symmetric quark-meson model was studied in Ref. (4| . 
The phase structures of the chiral quark-meson model and the Polyakov-quark-meson(-diquark ) m odel for 
the two lightest flavors as well as for 2 + 1 flavors were determined in Refs. T7, 48, 52. 53. 54. 55. 

Ellsl. ?7(l)A-violating terms in the FRC were introduced in Ref. [5|. For massless up and down quarks 
and a physical strange quark, the chiral transition was shown to be of first/second order without/with the 
U{\)a anomaly in the framework of a quark-meson model [HI, [6|. Also non-vanishing external magnetic 
fields were studied in such a model [Hllj El ■ Recent insights into second-order and fluctuation-induced 
first-order phase transitions in a C/(2)h x U{2)l symmetric model with scalar and pseudoscalar mesons 
were provided e^in Refs. [HHl E^- The influence of the U{\)a anomaly in a purely mesonic model was 


studied in Ref. [T|. The infrared-stable U{2)v x U{2)a fixed point detected in Ref. 0 was also found 
with the help of the FRC technique, but the subsequent stability analysis remained inconclusive [2(ij| 
FRC studies based on QCD degrees of freedom were presented in Refs. 65, H, 67, HI, 6| . 

Vector-meson fields were introduced into a chiral effective nucleon-meson Lagrangian in Refs. [zSInI 


lillTi El and the phase diagram at nonvanishing net-baryon density was studied within the FRC. 
However, the vector-meson degrees of freedom were only considered to be background fields, i.e., their 
fluctuations were neglected in the FRC flow. A first FRC study of (axial-)vector mesons in QCD at 
T = = 0 was performed in Ref. [7Hj |. revealing that the hadronic contributions to the flow in vacuum 

are dominated by the sigma meson and the pions. 

A study of the chiral transition at nonzero temperature within the FRC including vector and axial- 
vector mesons is, however, still missing. In this work, we fill this gap by applying the FRC to a purely 
mesonic model. We ascertain the mass deg enerac y of chiral partners above the transition and identify its 


order. We consider the so-called eLSM [Ty, [77|, [78| as effective model for the strong interaction at nonzero 


T and for vanishing /i. This model has been shown to re pro duce hadronic vacuum properties such as 
masses and decay widths to a surprising degree of accuracy [78l |. Furthermore, it has the same low-energy 
limit as QCD [79j. We study this model in a limit where it resembles the time-honored Sakurai model 

8|,E|. 


This paper is organized as follows: in Sec. 12.11 we introduce the two-flavor version of the eLSM. The 
subsequent part (Sec. 12.21) discusses the FRC formalism to describe fluctuations of spin-zero and spin- 
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one mesons. The numerical results for vanishing quark masses in the absence (or presence) of the axial 
anomaly are presented in Sec. 13.11 (13.21) . The case of explicit breaking of chiral symmetry by nonzero 
quark masses is discussed in Sec. 13.31 Section 2] concludes this work with a summary of our results and 
an outlook. 

We use natural units h = c = ks = ^ and work in a finite (3 + l)-dimensional Euclidean spacetime 
volume V x (0,1/T] at nonzero temperature T with periodic boundary conditions and, consequently, a 
discrete momentum spectrum: q = (a;„, q), where the Matsubara frequencies for bosonic fields are given 
by LOn = 2n'KT. We use a shorthand notation for spacetime integrations: 



with V = V/T. We employ Einstein’s summation convention, i.e., indices appearing twice are summed 
over. If these indices are Lorentz indices, ^ = 0,1, 2, 3, it does not matter whether they appear as co- or 
contravariant indices, because we work in Euclidean spacetime. We shall always use covariant Lorentz 
indices. 


2 Methods 


2.1 Extended Linear Sigma Model 


At low energies, quarks and gluons are confined inside hadrons, which are thus the effective degrees of 
freedom. An effective theory for hadrons must incorporate the chiral symmetry of QCD, as well as its 
spontaneous breaking. We work with a mesonic linear sigma model 82|, l83l | as an effective implementation 
of the strong interactiom which, besides scalar and pseudoscalar mesons, also includes vector and axial- 


vector mesons 


TdM 85 


the so-called eLSM [771 ■ The scalar and pseudoscalar fields are the 


real and imaginary parts of a complex Nf x Nf matrix E that lives in the [A^jl, Nf] representation of the 
group U{Nf){i X U{Nf)L- Under transformations of this group, E behaves as follows: 


U^EUl , 


( 2 ) 


where the group elements Ur^l are unitary matrices. In terms of hadronic fields, E = ((Tq -I- iTra)ta, 
with the generators ta of U{2) in the fundamental representation (tT[tatb] = Sab/‘2)- Here, Ua and tt^ 
represent scalar and pseudoscalar degrees of freedom, respectively. Analogously, we define right- and 
left-handed fields for (axial-)vector mesons (parametrized by axial-vector fields Aa^fi and vector fields 
Va,tj.)- Rfi = {Va,f, -f Aa,f,)ta, = (K,^ - Aa,f,)ta. They transform as: 


R/j. U\^R^UR , —>■ U\^LfjUR . 


The globally chirally symmetric Lagrangian is given by 


77|: 


£ = tr 
-btr 


( 7 ?^E)tz?^E] +mgtr(EtE)+Ai [tr + X^tr +-tv {L% + R%) 


+ A ) {Ll + Rl) - tr [H (E + E+)] - (det E + det E+) 


-ig 2 {tr{Lf,^ [L^,L^]} + tr{i?^,, + y H (E'I'E) tr -b i?^) 

-b/i 2 tr ^|i?^E| -b |EL^| ^ -b 2/13 tr (^EL^E^i?^) — [tr (L^LjyL^Lj/) -b tr {RfiRvR^Rv^\ 
[fi* “b tr {^R^R^RifR^y^ g^ tr (L^L^) tr {^RyRy'^ 

-ge [tr tr + tr (i?^i?^) tr (R^Riy)] , 


(3) 


(4) 


with the covariant derivative I?^E = 9^E -b igi{X2L^ — i?^E) and the right-/left-handed field-strength 
tensors — d^R^t and = d^L^, — d^L^. The term det E -b det E^ accounts for the U{1)a 

anomaly by breaking U{Nf)R x U{Nf)R to SU{Nf)v x SU{Nf)A x U{l)v- Its strength is determined by 
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the coupling constant ca- The flavor-diagonal terms tr \H (S -|- E1)] and tr [A (L^ -|- Rf^)] correspond 
to explicit symmetry breaking (ESB) in the (pseudo-)scalar and (axial-)vector sector, respectively: 


H 


Nf 

E^o 

2=1 


A = diag [(5 i,(52, ■ • ■,SNf] oc diag 




di ■ 


^Nf 


(5) 


For nonzero and degenerate quark masses = ... {Hq ^ 0, while all other Hq vanish; A oc 1 

has no further impact), these terms break the U{Nf)ji x U{Nf)L symmetry to U{Nf)v- 

It should be mentioned that there is a second way of introducing spin-one degrees of freedom to this 
effective theory. Within the gauged linear sigma model (gLSM) [8^ , (axial-)vector mesons are treated 
as massive Yang-Mills fields, accounting for the phenomenon of vector-meson dominance 8^, This 
model is constructed by requiring local U{Nf)fi x U{Nf)L symmetry. The gauge principle calls for a 
universal coupling of right- and left-handed vector fields to (pseudo-)scalars as well as among spin-one 
fields themselves. But due to the nonzero mass of the “gauge bosons”, of course, the gLSM is not a true 
gauge theory and the local invariance is already broken down to a global one. Since chiral symmetry 
is of global nature in QCD anyway, it seems to be logical to work with the Lagrangian ffl . Moreover, 
the “local” version does not reproduce the correct phenomenology of p and oi mesons [1, |83| , a problem 
which is solved by the globally symmetric eLSM [78|. 

In this study, we restrict ourselves to the isospin-symmetric two-flavor case, i.e., up and down quarks 
have the same mass. Hence we are dealing with scalar fields (cr, do), pseudoscalars (jj, tt), and with (/i, di) 
as well as (u;,;d) in the (axial-)vector mesonic sector. The fields cr, 77 , w, and /i are 5'f/(2)-singlet states, 
whereas the others form isospin triplets. The field p does not correspond to the physical p/p' mesons, 
which are mixtures of nn and ss (n = u,d stands for the nonstrange up and down quarks, s for the 
strange quark). As a first step in applying the FRG to the eLSM, we want to keep things simple and 
set the constants 52753 , ,56 as well as di, ^ 2 , da to zero. The above defined complex 2x2- matrices 

explicitly read: 


E = (cr -I- ip)to + (do + Of) ■ t , 

Rfj. = + {Pfj. + di^) • i 

T iPfi ai/x) ■ t 

and the different parts of C can be expressed as: 

'll 


( 6 ) 

(7) 

( 8 ) 


tr 


D^E = ^ [a^cr -b 51 {pfif, + TT ■ di^)]^ -b i [df,p - 51 (cr/i^ -b do • di^)]^ 

I 2 

+ 2 + di (Pm X do -b pSi^ + tt/i^)] 

1 2 

+ ~ 9i{'^ ^ PtJ. +craif, +aofif,)] , 


TOq tr (E'I'E) = ^ ^cr^ -b -b 5 ^ -b 7?^) , 
Ai [tr(EtE)]" = ^(a^+dB^ + 5^ + 7f2)' , 


A 2 tr 


(EtE)^ 






-b2 (cr^ -b 7f^) ( 5 ^ -b dp) -{ap-TT- dof I 


^ tr + Rl^) = ^ - d^uj^f -b ^ - d^p^f 

+ ^ - dyfifjf -b i {d^div - dydif^f , 

^t,{Ll + Rl)=^{fl + dl^+uTl + f^) , 

tr [id(E + Et)] =d°a , 

CA (det E -b det E^) = ^ ^cr^ — do^ — 5^ + • 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 
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Obviously, the w meson completely decouples from any interactions. This would not be the case had we 
included other terms from Eq. 

In the low-temperature broken phase, the isoscalar a field acquires a non-vanishing expectation value 
(tT)o = {qq) = const. ^ 0 (here, the angular brackets are the notation for expectation values and the 
subscript 0 denotes the absence of external sources). Therefore, one has to consider fluctuations around 
the physical ground state and thus a shift of the a field by its expecation value: a ((7)q -I- a. The 
expectation value (fT)^ acts as the order parameter for the chiral phase transition. After accounting for 
this shift, integration by parts then gives rise to the bilinear terms gi (a}g and gi {a)^ tt ■ d^di^. 

They represent the so-called tt-oi- and 77-/i-mixing, leadin g to nondiagonal elements in the scattering 
matrix. Usually these terms are eliminated. Following Ref. [SlJ, this is done by shifting the axial-vector 
fields: fi^ -)■ /i^ -|- wd^i-q and dif^ -)> -I- with w = gi (ct)^ /[m\ + {gi (cr)o)^]- In turn, the 

axial-vector fields become explicitly RG-scale dependent (through the dependence of w on (cr);,) and the 
pseudoscalar states need to be renormalized: —>■ v^^TTa, Zt^ = l-\- (gi {(j)^)^/m\. This provides the 
canonical normalization of all one-meson states, such that their Fourier components can be interpreted 
as creation and annihilation operators in the process of quantization . For a precise discussion of the 
a shift and its implications on the FRG flow we refer to Ref. [t^. Instead of redefining the oi and fi 
fields, one may also work with nondiagonal propagators as performed in Ref. [s^ • 

The vacuum expectation value ((T)g is the minimum of the classical potential energy density U((cr)): 


^ (^0 - ca ) {erf + ^ (^^1 + y ) ~ ’ 


dU 


d(CT) 


= 0 . 


{<r) = {<r)o 


(17) 


The wave-function renormalization Z.^ is related to ((T)g and to the masses of the ai and p mesons by: 


(cr)o = 


TUa 


(18) 


where /tt — 93 MeV is the pion decay constant. The KSFR (Kawarabayashi-Suzuki-Fayyazuddin-Riazud- 
din) relation 88, 8^ predicts that Z^^ = 2. This slightly differs from the value of ~ {rna^jmpf ~ 2.552 
quoted by the particle data group [^ . 

The identification of the mesonic Helds with measured resonances listed in Ref. is partly straight¬ 
forward: the pions and the q (as the pure fin state arising from unmixing the physical q and q') have 
a mass around 140 MeV and 700 MeV, respectively. The vector fields w and p represent the a;(782) 
and p(770) resonances. The axial vectors fi and ai correspond to the /i(1285) and ai(1260). For 
the a and the gq fields, however, it is controversial whether they should describe {/o(500),ao(980)} or 
{/o(1370), ao(1450)}. It was argued in Refs. that the latter option might be favored. 


2.2 Functional renormalization gronp 

The Wilsonian renormalization group (RG) performs the mode integration of (quantum-)statistical fluc¬ 
tuations from the ultraviolet (UV) to the infrared (IR) in a stepwise manner, i.e., it successively takes 


momentum-shell by momentum-shell into account |9lL l92 


. a stepwise manner, i.e., it successively taxes 
H [11 m,S111- The FRG is an implemen¬ 


tation of this procedure which allows us to non-perturbatively formulate quantum field theories in terms 
of a differential equation. This flow equation dictates the scale (k-) dependence of the effective average 
action T/j, which interpolates between the bare interactions at some UV cutoff scale kjjv = A and the 
macroscopic physics including all fluctuations in the IR, /cir = 0. A /c-dependent term ASk is added 
to the action S in order to provide an effective cutoff at momenta q^ ~ k^, such that only modes with 
q^ > are integrated out in the RG flow. The term ASk regulates the scale evolution of Tfc in such 
a way that the full effective action T = Tfe-^o is obtained in the IR limit. The effective action T is the 
generating functional of one-particle irreducible vertex diagrams of the theory. For fc —>■ A, in contrast, 
the classical action is recovered: Tfc^A = S. 

Our investigations focus on the theory defined by the Lagrangian ([4]). In order to simplify the following 
discussion, we denote (pseudo-)scalar fields by (pi and (axial-)vector fields by Ai^^. The field-strength 
tensors of the latter are denoted as = dpAi^^ — d^Ai ^. The fields pi and are subject to thermal 
(for r > 0) as well as to quantum fluctuations (also at T = 0). Following the discussion in Ref. 98|, we 
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apply Stueckelberg’s Lagrangian 9^ lOPjl with coupling Ast to derive the FRG flow equation: 


S= C 




( 19 ) 


(Axial-)vector mesonic fields usually have three physical degrees of freedom. The additional term in 
Eq. m, however, promotes the unphysical fourth to a physical one. Hence, in the following, all vector 
helds initially have four instead of three degrees of freedom. Although not necessary to ensure renormal- 
izability (^ . this formalism guarantees that we work with invertible inverse tree-level propagators [sij. 
Furthermore, this strategy allows to derive the grand canonical partition function in a simple manner 

M- 

As a starting point for deriving the FRG flow equation of the theory at hand, we consider the scale- 
dependent generating functional H4 for connected Green’s functions: 


Wk [Ji, Jz,k.] = In Zk [Ji, J 




( 20 ) 


As discussed above, we add a regulator term lS.Sk to the action: 

which can be interpreted as a momentum-dependent mass term, and we also included sources Ji and 
for scalar and vector fields. To ensure the required UV/IR limits for the flow of F^, the regulator 
functions R%{q) and R^ must fulfill the following relations: R%{q)^ R^f^^{q) —>■ 0 for fc —>• 0 , as well as 
—>■ oo for fc —>■ A. On top of that, the regulators should satisfy: i?®((z), ^^^ 1 ,( 9 ) for 
q 0 and .Rfc(9), ^^(?) ~ 0 for g — 00 . Apparently, low-energy fluctuations {q^ k^) are effectively 

separated from the RG integration process by giving them an additional “mass” ~ while fast modes 
(g^ ^ fc^) are not influenced. The effective average action F^ -|- ASk is the Legendre transform of Wk, or 
in other words: 


Ffc [4^k,i, Ak.i.k.] — E ^ ) [Jii q)^k.i{si) A ?)^fc,2,^(?)] 

9 


A^k [Ji, ASk [4^k,i, Ak.i.fi] , 


( 22 ) 


where (^fc,i(g) = (‘/?i(g)) = V ^6Wk/SJi{-q) and Ak,i,f,iq) = (^i./x(?)) = ^ ^SWk/5Ji^f,{-q) are the 
expectation values of the fields in the presence of the sources Ji and Ji^^. Although the Legendre 
transform Ffc -|- ASk is convex by definition, this does not hold for Ffc itself, as ASk is not necessarily 
curved in the same way. Exclusively in the case fc —>■ 0, where ASk —^ 0, Ffe-^o = F becomes the true 
Legendre transform of H4_>o = W, and thus is definitely convex. 

For fixed values of the fields, differentiation of Eq. (12^ with respect to k yields the FRG flow equation: 

dkTk = ^ {tr [GS(g, q)dkRtiq)] + tr [GV^,(g, q)dkRl,,^iq)] } . (23) 

q 


Here, G® and G^^,^ denote the full propagators for scalar and vector fields. Introducing the general field 
notation $ = {(pk.i, Ak^i.n) and using the fact that VG^ = -|- , Eq. (1^ simplifies to: 


SfcFfe = 2 fo 


(rf+Rfc) 


-1 


^k'R'k 


(24) 


The momentum summation has been included in the definition of the trace. The propagators Gfc and 
Rfc are matrix-valued in momentum space and in all internal spaces. 

In principle, the FRG equation and the resulting macroscopic physical observables should be indepen¬ 
dent of the form of the regulators, which are only restricted by the limits discussed above. In this case, all 
trajectories in coupling space predicted by different choices of Rfc start at the point rfc_i.A = S and ter¬ 
minate at Ffc_).o = F. In practice, however, one needs to truncate the infinite hierarchy of flow equations 
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arising from Eq. (l24l) in order to solve them. Indeed, this fact inevitably leads to a regulat or-d ependent 


bias, which, fortunately, can be minimized by working with the optimized Litim regulator 101 


One convenient truncation scheme is the expansion of Tk in terms of field derivatives [^, l43| . For a 
purely scalar theory, this would read: 


rfc[<&] = [ 

J X 


t/fc($) + + o{d^) 


(25) 


C/fe is the scale-dependent effective potential, symbolizes the wave-function renormalization associated 
to the scaling of the kinetic part. The LPA assumes the wave-function renormalization in the derivative 
expansion (1251) to be field-independent and fixed to its initial value of one, Zk{^) = Zfc = 1, which 
corresponds to a vanishing anomalous dimension. Momentum-dependent interactions are neglected in 
LPA. In the case of the eLSM, the second term in Eq. (ESI) is given by the kinetic terms of spin-zero and 
spin-one fields (modified by Stueckelberg’s Lagrangian): 


= f 

J X 






(26) 


— d^Ai^fj,. At nonzero temperature, it is technically adva ntageous t o employ the three- 
dimensional version of Litim’s optimal regulator oc (fc^ — (f^)0(fc^ — (f^) [l02l [l03|. Remembering the 
above discussion about their various limits, the regulating functions are chosen as: 

Rl(q) = {e - g2)0(fc2 _ ^-*2) ^ (27) 

= [nj.(g) + Astn;(,(g)] {e - - f) , (28) 

where we have defined the transversal and longitudinal projection operators nj,^(g) = — q^q^/q^ 

and n^j,(g) = q^q^/q^, respectively. Since dk^k = VdkUk for spatially uniform field configurations 
^i{x) = and since in the limit E —>■ oo we have l/EX^q I d^g/(27r)^, using Eq. as well as 
Eq. (1^ . Eq. (1^ becomes: 


dkUk = TY^ 


f y r 

Jvik){2nf \ 


1 trG®(a;„ 


. 9 ) + [nj>.(9) 




(9)] tr 


(Wn,g)| , 


(29) 


where G = VG and V{k) denotes the spherical volume with radius k in three-momentum space. 


3 Results 


In this section, we numerically solve the ERG flow equations in LPA for the eLSM introduced in Sec. 12.11 
for three scenarios: (i) the U{2)fi x f/(2)i-symmetric case, (ii) the case with U{1)a anomaly, and (hi) the 
case with U{1)a anomaly and ESB. To this end, we discuss the expansion of the potential Uk in terms of 
the respective invariants under the given symmetry and fix the bare couplings in the UV such that the 
renormalization flow produces reasonable values for the physical observables in the IR. From the behavior 
of the chiral order parameter ((t)q as a function of temperature we infer the order of the phase transition 
and illustrate the restoration of chiral symmetry by computing various mesonic screening masses. 

Let us remark that there are, in principle, two different strategies to proceed [^. In the first strategy 
one expands the potential Uk around a (local) minimum. The advantage of this method is that one has 
to solve only a few flow equations (one for each coupling and an additional one for the scale-dependent 
order parameter). In doing so, however, we can only deduce the potential right at a local minimum. It 
is not clear whether this local minimum is also the global one. Especially in the case of a first-order 
transition with two emerging minima, it is crucial not just to know the potential at the expansion point 
but also at any other local extremum. 

To overcome this difficulty, in this paper we follow the second strategy, where Uk is discretized on 
a grid. Here we gain information about the entire form of the potential, but this strategy needs a 
lot of computational power as we have to solve flow equations for each grid point. Nevertheless, we 
utilize this approach because it allows us to figure out the transition order in a comparatively quick and 
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uncomplicated fashion. We tune the potential in such a way that the physical configuration is located at 
(cr)p = v ^/77 for fc 0 . 

Remember also that the effective action Ffc is a functional of the classical fields = {<fi) and 
of. Eq. (1^ . but for the sake of simplicity the brackets indicating expectation values will 

be omitted in the following, e.g. {rj) —>■ r]. Another point is that we are setting certain fields to zero after 
( 2 ) 

the calculation of , since we only need to consider as many fields to be non-vanishing as there are 
independent invariants. 


3.1 Chiral limit without anomaly 

For zero quark masses and in the absence of 17(l)A-symmetry breaking (/iq = 0 as well as ca = 0), 
the Lagrangian (U) is invariant under the full U{2)ii x U{2)l symmetry. In the LPA neither wave- 
function renormalization nor momentum-dependent interactions are taken into account. Since the in¬ 
volved (axial-)vector mesons typically have a mass of < 1.2 GeV, an ultraviolet cutoff of A = 1.2 GeV is 
chosen. After substituting the fields by their expectation values, in compliance with Eq. (l26)) we derive 
from Eq. the following effective average action at the UV scale: 

11 11 1 2 I 2 

-5^cr9^cr -I- -d^do ■ dfj,do -I- -df_tr]d^r] -|- ■ 9^7? -|- - {dfj,uj„ - d^ujfj.) -I- - - d^Pfj.) 

+ ^ -t- i {dfj,di^ - d^dif^f + ^ 

+ Ci,aC + C2,aCi a C3,A^2 + C4^A^3 + C5^aC4 C • (30) 





In Eq. (1301) we have introduced the 0(8) mass invariants C s-nd ^4 well as the other 17(2)/j x U{2)l- 
symmetric expressions ^ 2 , ^3 contained in Eq. ®. They are linear combinations of the different interaction 
terms in the effective potential, namely: 


= a'^ + dl + T]'^ + , 

6 = {<y^ + 7 ?^) + Qq) - (o-p - 7? • oo)^ , 

^3 = (7? • dif, + -F (oo • oi;, + cr/i/i)^ 

+ (pfj, X dQ + pdifj. + -I- (tt X -I- adifj. + oo/i/i)^ 

_ n 2 I ~*‘2 I 2 I ~Q 

?4 — fifj, + +^tj.+ Pfj. ■ 

The scale-dependent couplings are defined as: 

k 

Cl,fc = , C2,fe = 


T 2'^2,fc 


C3,fc = 


^2,fc 


C4,fc = 


alk 





(31) 

(32) 

(33) 

(34) 


(35) 


For the expansion of Uk we have to replace all non-vanishing field variables by appropriate expressions 
of the invariants {Ci^ 2 )Csi f 4 }- In order to do this, we keep the fields a, aj, oJq, and pj nonzero. All 
others are set to zero. Solving the four Eqs. (ED) - (El for these four variables yields: 
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2 



1 ({1 + vVbfl) 

§3 

c ’ 


(“0) “ 2 

W)' = e4-|. 


(36) 


At first glance, this mapping seems to be singular for —» 0 , b ut we have checked that all parts of 
the flow equations oc cancel for = 0, cf. also Refs. 6 ^, 104 1. Furthermore, the mapping does not 
preserve Euclidean invariance, since we keep only the p = 0-component of the vector fields p^ and a\^. 
This gives rise to unequal screening masses of some vector components with differences oc ^ 2 , ^ 3 , or ^ 4 , 
but this is also not relevant since we assume that only the sigma field acquires a non-vanishing vacuum 
expectation value (and hence —?► ^10 = ctq and ^ 2 , ^ 3 , ^4 —t 0). This means that we are concerned with 










a one-dimensional investigation along the ^i-axis and that, in this limit, Euclidean symmetry is restored 
again. 

Without axial anomaly and quark masses, we expect the chiral phase transition to be of first o rder as 
argued in Ref. 0 (among many other studies) and summarized by the “Columbia plot” 


105, lOa. Since 


^ 2 ) and ^4 are set to zero in the end, it is reasonable to truncate Uk at linear order in these invariants 
(although the flow equation generates terms of arbitrary order in these invariants), with coefficients that 
are functions of 


C^fc(6j?2,C3,'C4) = 14 (^ 1 ) -I-Wfc(^l)^2 +-^fe(Cl)C3 + ^fe(6)?4 ■ (37) 


The physical vacuum is specified by the condition dUk/da = 0 for a = ao, and the squared mass of 
the a field is identical to the curvature of the effective potential: 


dUk 


da 

d^Uk 


da^ 


= 2v^^fe(6o) = 0 , 

= 2V'fe(6o) + 4fiol4"(.^io) ■ 


(38) 

(39) 


From Eq. (1381) one sees that, for ^10 ^ 0, the minimum of the effective potential can also be determined 
from the condition = 0. Once the system changes to the restored phase (ctq = ^10 = 0), can be 
nonzero, though. 

The explicit flow equations for I 4 , Xk, and Yk are obtained by differentiating Eq. (l29)) and 
evaluating it for the physical configuration: 


dkVki^i) — ^fet7fe(^i,C2,C3j'C4)l52,{3,{4=o — 2 ’ 

dkWki^i) = 9fe%C/fc(6,e2,6,e4)l5,.53,5,^0 ^ 

dkXki^i) = 5fe%t/fe(ei,6,6,e4)|j3,53,^,^o = (-Gferg3GfeafcRfc 

dkYki^i) = 5fe%[/fe(6,e2,6,^4)153.53,5,^0 = (-Gferg,GfeafcRfc) 


52.^3,44 = 0 
52 . 53,54 = 0 
52 . 53,54 = 0 


(40) 

(41) 

(42) 

(43) 


The equation for Vk turns out to be equivalent to the flow equation for free fields, as interactions are no 
longer present for 6,6,6 = 0: 


Tk^ 


atUKi)=g^i: 


1 


u;2 -h fc2 + 2Vl wl + k‘^ + 2Vl + 46 14 " ujl + k‘^ + 2-b 26 Wk 

12 12 
-b 


a;2 -b fc2 -(- 26s ^2 -b A:2 -b 2(6; + 6^fe) 

4 4 


a;2 -b fc2 -b 26/Ast ^2 -b A:2 -b 2(6 + 6^fe)/Ast 


(44) 


Here, explicit 6-dependences on the right-hand side are implied. The masses of the fields p and w as well 
as /i and oi are equal in the presence of SU{2)v invariance. Equation (Hll) reveals some characteristics 
of the considered model: the 4 -b 1 -b 3 = 8 (pseudo-)scalar degrees of freedom are represented by the 
first line. The first term of Eq. corresponds to the mass-degenerate p and if. The second and the 
third describe the a and the do, respectively. As a cross-check, note that the mass eigenvalues of these 
spin-zero fields are equal to the ones quoted in Ref. [^. The second line corresponds to the physical 
degrees of freedom of the (axial-)vector mesons ((4 x 3) -b (4 x 3) = 12 -b 12 fields), whereas the last 
line corresponds to their unphysical degrees of freedom (4 -b 4) introduced via Stueckelberg’s Lagrangian. 
Obviously, these eight additional degrees of freedom decouple from the flow for Ast —>■ 0. The Mats ubara 
sum over n can be carried out analytically, e.g. by a contour integral in the complex plane |l07l |. The 
flow of 6(6) entangles with the flow of VEfc(6), Aife(6), and 6(6)- Thus the flow equations for these 
coefficients are necessary to obtain a closed set of differential equations. For the sake of clarity, they are 
presented in Appendix lAl 
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Figure [T] A summarizes the dependence of (Tq on T. The order parameter (Tq becomes successively 
smaller and drops discontinuously to zero at Tc ~ 147.4 MeV, indicating a first-order phase transition. 
The discontinuity occurring at this temperature is marked with a dashed line. Figure [T]B demonstrates 
how, as the temperature increases, the mesonic screening masses of chiral partners approach each other. 
These are (i) p and oi, (ii) uj and fi (their masses are not shown explicitly, since = rup and = 
ma^), (hi) a and tt , as well as (iv) oq and 77 {rrip = 171-^, thus we do not show nirj explicitly). The masses 
of chiral partners become degenerate at the transition temperature and above. Note that solid lines 
correspond to data interpolated using cubic splines. In cases where the data points are explicitly given in 
terms of colored crosses, however, the lines represent a cubic smoothing spline ht. Details are provided 
in Appendix IbI 

The p/oi mass increases/decreases before the transition point is reached. Pions and ry are the Goldstone 
bosons of chiral symmetry breaking, and thus necessarily massless in the broken phase. This is explained 
by Eq. (1551) : if ^10 ^ 0, the physical minimum is located at the point where V)((^i) = 0. Inspecting 
Eq. ((551) . we see that I 4 ' is indeed proportional to the (squared) mass of pion and 77 . However, for T > Tc 
the ground state is characterized by ^10 = 0 , and thus I 4 '( 0 ) (i.e., the masses of pions and 77 ) may differ 
from zero. 

We tuned the UV-parameters in the vacuum to achieve the most “realistic” mesonic screening masses 
and a nonzero value for (Tq of around 147.9 MeV (the tree-level value is tro = /tt — 148.8 MeV). Ap¬ 
parently, the IR vacuum masses of the a and oq mesons are far too small compared to what we expect from 
the PDG [ 9 ^, no matter whether we choose the assignment {/o(500), ao(980)} or {/o(1370), ao(1450)}. 
We will return to a discussion of this issue below. Furthermore, the masses of w and p are too heavy. 
The ratio between the masses of p and oi is smaller than expected (experimentally it should be around 
1 . 6 ). 


A B 




Figure 1 : Phase transition in the eLSM with t/( 2 )_R x U{2)l symmetry. A Order parameter ao as a function of 
temperature. B Screening masses as a function of temperature. 


3.2 Chiral limit with U{1)a anomaly 


As a second scenario we want to study how the U{1)a anomaly influences the order of the transition 
and the mesonic masses. The coupling ca in Eq. (|5]) is now nonzero and quantifies the strength of the 
17(l)A-symmetry breaking. Proceeding similarly as above. Fa is slightly modified: 


rA = 


J I ■ dpQo + ^9,^7719^77 + ■ dpTT + i {dpUJi, - duLOpf + i {dpPi, - d^ppf 

+ ^ {dpfii, - d^fipf -f i {dpai„ - d^aipf + ^ {dpWpf + {dpPpf + {dpfipf + {dpdipf 


+ Cl.A'Cl + C2,k£a + C3,aC'?2 + Ci,K^2 + Cs^aCI + C6.A'?3 + C4,aC3 + C5.aC4 ( , 


(45) 
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with the new invariants ^ 2 , ^ 3 : 

(46) 

(47) 

(48) 


^1 = cr^ + 7?^ , 

^2=rf + al , 

6 = (o’^ - 7? • a^f . 


The former invariants and ^2 are functions of the new invariants ^ 1 , ^ 2 , and ^ 3 : + ^ 2 , ^2 = 

^ 1^2 ~ ^ 3 - The origin of the new invariants is the J7(l)^-symmetry breaking term ~ det S + det in 
Eq. (|4|). The invariants ^3 and ^4 remain unchanged. The scale-dependent couplings of the effective 
potential are now defined as follows: 


ClM = 2 
1 / 

C4,fe = 2 V™ 

_ 5i.fe 

C4.fe- 2 


0 ,fe “ C.A,k) , 

C2,/c = 

i' 

^Ai^fe -b —A2,fc^ 

5 

O.fe + CA.fc) , 

C 5 ,A: — 


( Ai^fe -b — A2,fc j 

•) 


C5,fe = 


1 

2 


^1,- 






(49) 


We again want to expand Uk in terms of ^ 1 ,..., ^ 4 . For five different invariants we have to keep at least 
five fields nonzero. In addition to cr, aj, oJq, and pj we decided to take the 77 field into account and map 
those variables onto the set C 2 ,Csi^ 3 ;^ 4 }- 


a 


2 


6 1 




6 
6 ’ 


6 ’ 



6 

6+6 


{plf = 


6 

6 + 6 


(50) 


Similarly to Sec. 13.11 the potential is now discretized with respect to and one assumes that the physical 
minimum is attained for ^10 = fig and 6,6 = 0. Thus, in the ansatz for the effective potential, we restrict 
ourselves to the linear order in ^ 2 , 6 , 6 , ^ 4 : 


c^fc(6,6,6,6,6) = Efc(,^i) -I- iTfc(6)6 + d^fc(6)6 + +^(6)6 + ^fe(6)6 ■ (5i) 


The flow equation for 14(^i) reads: 


'Thi 

aM&) = 


+ 


1 


[ujl + k^ + 2Vl ujl + k^ + 2Vl + 4^iVl' 
1 12 


-b 


-b 


ujl + k'^ + 2Wk 
12 


-b 


+ ^^ + 2(W4 -b 6 "^fc) -b -b 2Zk uj'^ + k"^ + 2{Zk -b 6 T/c) 
4 4 


-b 


+ ^^ + ‘2‘Zk /Ast -b -b 2{Zk -b 6 Tfc)/Ast. 


(52) 


Note that the rj meson now attains a mass different from that of the pions: it becomes massive in the 
spontaneously broken phase, as it should be, since U{\)a is explicitly broken and there are only three 
(and no longer four) Goldstone bosons. 

Again we tune the bare parameters such that most realistic meson masses are obtained in the IR. The 
order parameter (Tq continuously decreases with increasing T, until it vanishes at a critical temperature 
of approximately 276.5 MeV, see Fig. [5] A. We conclude that the axial anomaly turns the first-order 
transition, found in the U{2)fi x 17(2)i-symmetric theory, into a second-order phase transition with a 
significantly higher critical temperature. Figure [2] B shows the evolution of the masses of (pseudo-)- 
scalar and (axial-)vector mesons. Again, there are four different pairs of chiral partners, (p,ai), (w,/i) (the 
masses of which are identical to the corresponding mesons in the first pair and thus not shown explicitly), 
(tT,7r), as well as ( 77 , 09 ). The masses of chiral partners become degenerate at the transition temperature. 
The pions assume a non-vanishing mass above the critical temperature for the same reason as discussed 
in the previous section. At zero temperature, the p{ui) meson mass is close to its physical value, but the 
ai(/i) mass is too small. The mass difference between vector (w, p) and axial-vector mesons (/i, oi) 
comes out to be ~ 268.7 MeV and is of the same magnitude as for the U{2)b. x {7(2)i-symmetric case 
(~ 265.4 MeV). The 77 -meson mass is too large, but remember that the 77 in our case consists only of 
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Figure 2: Phase transition in the eLSM with SU{2)v x SU{2 )a x U{1)v symmetry. A Order parameter cro as 
a function of the temperature T. B Screening masses as a function of temperature. 


nonstrange quarks, while the physical 7] is an admixture of nonstrange and strange quarks. The vacuum 
mass of the a of around 357.4 MeV is now only slightly smaller than the experimental value for the mass 
of /o(500), while the vacuum mass of uq is very close to its experimental value. The data points for the 
mass of the a fluctuate strongly as a function of temperature, since the potential is rather flat in this case 
and its curvature (the squared a mass) is rather hard to determine numerically with reasonable accuracy. 

3.3 Explicitly broken chiral symmetry with anomaly 

Finally, we apply our FRG analysis to the case of ESB due to nonzero and degenerate quark masses 
{Hq ^ 0). The truncation of the effective potential (IFTl) is still valid. In the case of ESB, the root of Vl 
no longer coincides with the one of dUk/da for nonzero vacuum expectation values of the cr field. The 
global minimum ^lo must now fulfill the following relation: 

2v^^fc(eio)-/i!S = 0, (53) 

i.e., the expansion coefficients Wk-, Xk, etc., are evaluated for a shifted fio, producing massive pseudo- 
Goldstone bosons tt (which have a nonzero mass even for uo 0). 

According to Eq. (I53L the potential never has a global minimum at = 0. Solving the flow equations 
on the grid, the results are shown in Fig. [3] In Fig.|3]A, we see that (Tq decreases and tends asymptoti¬ 
cally towards zero. Gonsequently, the transition is a crossover transition. At an estimated pseudocritical 
temperature of Tpc ~ 354.8 MeV, the curvature changes its sign. Figure[3]B reveals that, with the choice 
hjj = 3 X 10® MeV®, the pions exhibit a nonzero mass of ~ 142.5 MeV in vacuum. The masses of chiral 
partners approach each other, but do not become identical. We further observe a dropping oi meson 
mass but an increasing p mass. The vacuum cr mass of 598.5 MeV is a little larger than the physical 
value for the /o(500) resonance. The remaining masses are in good agreement with the results of Sec. 
o The gap between the masses of the p and oi mesons is now 314.5 MeV (compared to 268.7 MeV in 
the previous case). 

4 Summary and Outlook 

The QGD transition is commonly associated with the restoration of chiral symmetry. Experimentally, 
this could be detected by a change of the in-medium masses of (axial-)vector mesons. It is therefore 
essential to include them in a theoretical analysis. Non-perturbative continuum methods, such as the 
FRG, provide new insights into the QGD transition, as they do not rely on weak couplings and are 
applicable at nonzero net-baryon density where lattice QGD suffers from the fermion-sign problem. 
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Figure 3: Phase transition in the eLSM with ESB and U{1)a anomaly. A Order parameter ao as a function of 
temperature. B Screening masses as a function of temperature. 


In this study, we investigated the chiral transition for two flavors by applying the FRG formalism to the 
eLSM, an effective low-energy model for QCD. Thus, our work is an extension of many studies involving 


two-flavor effective models for QCD, see e.g. Refs. [18|, l63|, [6^ , in the sense that vector and axial-vector 


mesonic degrees of freedom are now incorporated into the FRG flow. In order to derive the FRG flow 
equations with (axial-)vector mesons, Stueckelberg’s Lagrangian has been employed 98, |^, IiO(TI| . We 
use the grid method and the LPA to compute the flow of the effective potential. The order of the phase 
transition and the meson screening masses were determined in three different scenarios: (i) the chiral 
limit without U{1)a anomaly, (ii) the chiral limit with U{1)a anomaly, and (iii) the realistic case with 
nonvanishing quark masses and U{1)a anomaly. 

Overall, our numerical results are broadly consistent with previous findings. Regarding the full U{2)rX 
C/(2)i-symmetric theory, cf. Sec. 13.11 our conclusion is compatible with Ref. [6J| and the statement of 
Pisarski and Wilczek ITj that the chiral phase transition is of first order for Nf = 2 and massless quarks. 
Reducing the symmetry to SU{2)v x SU{2)a x U{1)v ISec. 13.21) turns it into a second- order transition, 
which is also in agreement with Ref. 0 as well as with the three-flavor study of Ref. [^. Explicitly 
breaking chiral symmetry to an exact isospin symmetry generates a crossover transition (Sec. I 

In comparison with the results from the CJT formalism or lattice-QCD simulations 
the pseudocritical temperature Tpc ~ 354.8 MeV of the crossover transition comes out larger, cf. Tpc ~ 195 
MeV resp. 155 MeV in the aforementioned approaches. It is conceivable that this deviation arises from 
the lack of quark fields in our approach, which, being rather light degrees of freedom, evidently contribute 
substantially to the FRG flow [75|, and from the fact that presently our approach ignores momentum- 
dependent vertices and interactions among (axial-) vector mesons. It has to be clarified how higher orders 
in the derivative expansion affect the transition. We hope that by additional investigations, as stated 
below, we are able to improve our results. 

In all three scenarios studied here, it was demonstrated how the masses of chiral partners become 
degenerate at the phase boundary and beyond, see Figs. B and[2]B. The mass degeneracy is a 

necessary condition for the restoration of chiral symmetry [y. III, ISl, ISj . Let us note that in our study the 
mass of the ai mass decreases towards the chiral transition, but not the mass of the p, cf. Fig. [3]B. In 
the GJT study of Ref. ^ the authors also found an increasing p mass towards the chiral transition. In 
principle. Ref. argues that the p meson mass has to increase in the framework of a gauged two-flavor 
LSM, but a globally symmetric LSM could also allow for a dropping p mass. 

In the physically most realistic scenario with ESB and U{\)a anomaly, the vacuum masses of (cr, oq) 
come out to be (598.5, 996.3) MeV. This is in the range of the masses of the light scalar resonances 


{/o(500),ao(980)}. However, Refs. [76|, [77|, [78| suggested that the chiral partners of tt and p should 


be {/o(1370), ao(1450)}. This seems to be a natural scenario, if the latter are (predominantly) quark- 
antiquark states. Then, the light resonances {/o(500), ao(980)} are most likely made of four quarks, e.g. 
in the form of resonances in the scattering continuum, or even bound states, of two pseudoscalar mesons. 
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Also bound states of diquark and an ti-di quark molecules have been suggested to explain their nature. As 
suggested a long time ago by Jaffe [l09l| . this ’’tetraquark” interpretation of the light scalar resonances 
would naturally explain the ’’inverse mass ordering” of these states. By construction, the FRG approach 
resums correlations of infinite order and thus, if {/o(500), ao(980)} are correlated states of pseudoscalar 
mesons, would naturally generate these mesons dynamically. This could be an explanation why the 
masses of a and oq come out close to those of {/o(500), ao(980)}. (In fact, we were not able to find 
UV-parameters such that the IR vacuum masses of cr and oq are close to those of {/o(1370), ao(1450)}.) 
Note tha t the chiral transition was studied in the presence of both a light and a heavy scalar state in 
Ref. [llOl| . There it was shown that there is actually no conflict with the ’’tetraquark” scalar state being 
light and the heavy scalar state being the chiral partner of the pion. 

There are many questions left open for future study. E.g. one should investigate the order of the phase 
transition as a function of the anomaly strength. The first-order transition in Sec. 13.11 should smoothly 
pass into one of second order, as shown in Sec. 13.21 Moreover, one needs to hgure out when exactly the 
U{1)a anomaly disappears for high temperatures. This can be done by assuming ca to be proportional 
to an explicitly T-dependent instanton density. In order to decide whether the transition lies in the 0(4)- 
universality class or not [HEl, the critical exponents have to be calculated. A natural next step in our 
analysis is to account for non-trivial wave-function renormalization factors, i.e., going beyond the LPA. 
One can readily extend our investigations to Nf = 3 quark flavors, but in this case one has an additional 
order parameter (the strange condensate) which necessitates the use of a two-dimensional grid (b^ and 
thus considerably increases the numerical effort. Baryonic degrees of freedom also play an important 
role for dilepton production M- Therefore, it is mandatory to involve them in the FRG flow. The first 


111 . Ill2l| - if the latter 


candid ate f or such an extension would be the nucleon and its chiral partner 

exists |ll3l| . As this partner would probably be hea vier than the A baryons, one should furthermore 
include spin-3/2 resonances and their chiral partners 114| . Finally, another topic for future studies is to 
introduce quarks in the FRG. The reason for this is that, as mentioned above, a large contribution to the 
FRG flow usually comes from the quarks, while the vector and axial-vector mesons are suppressed due 
to their large mass. 
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A Flow equations 


The flow equations for the expansion coefficients of the effective potential Uk are generally of the form: 


^kFk,^ = ^ E / fhAFk,FiFl!,Lol, f) 


(54) 


with Fk = (14, IFfe,...), fk = {vk,Wk, ■ ■ ■), and q = |g|. Dependences on were omitted. The 

derivatives = dFk/dA and F" = d^Fk/id^? (or F^ = dFk/dA and F" = d^Fk/idA)"^, respectively) 
are approximated by finite differences, which leads to a solvable system: 

dkFk,i = Q'^fk,i{Fk,ujl,q'^) . ( 55 ) 

n ^ 

The functions fk,i are specified in subsections lA.ll and IA.21 The sum over Matsubara frequencies is 
performed analytically, e.g. 



1 

F^ + qA 



coth (^) 

2gF(fc2+4aK' + 2R,'-<72)2 


14 









erf 

■ 8r2 (fc2 + 4^^^" + 2V^; - g2): 


g^csch 


+ 


2 ^ Vfc2+46V;;'+2V^' 


2T 


+ 


8r2 (P + 4^^^. + 2^,') (fc2 + 4^1 y" + 2V' - q^) 2 
8T (i;2 + 4{,14'' + 2V2) V" (t= + 4{iVJ' + 21^ - ,2) 2 

3,i„,. UE±^!^ns\ <,,d,2 


STy'k^ + 46 + 2Vl (F + 4^1 y" + 214 ' - 9") ■ 


(56) 


Where necessary, the momentum integration is performed by using numerical quadrature. 

The potentials 14 and 14 are initialized as follows: 

h4(fi)=aA(ei-5i)^ (57) 

Mii) = -aK{h-bl)\ (58) 

where oa = 4.0, 6a = 297.6 MeV, Sa = 5.5, and 6 a = 255.0 MeV. Furthermore, we have 114v(Ci) = 30, 
Xa(6) = 20, 14(6) = 2.87 X 10^ MeV^ Wa( 6) = 2.5 x 10^ MeV^ Xa( 6) = -1-0, 14(6) =31.0, 
and Za = 1.27 x 10® MeV^. Under the assumption of ESB, the UV parameters change to: 1Fa(6) = 
2.0 X 10® MeV^, 14(6) = 45.0, and Zp^ = 1.1 x 10® MeV^, while the rest remain identical. 


A.l Flow equations without U{1)a anomaly 

In the case without U{1)a anomaly, the quantities fk,i in Eq. (IKKl) read 


, / 3 12 12 1 4 \ 

if 4[214"(8Wfe-361U^) + 156W'fclU'+46^1^^+ 1^2] 2AXk 192^2 

= n-+ ^ + ^ 

4 [2U4 {2e,Wi: + 96+ 61Ufc) + 61^^ (31U^ - 26VF;) + 46lU^2 ^ 3p^,2] 

6E4 (2U4 - W,) 

8(261U^ + lUfc) 161^2 24X4 

E® 6£^p^ J ’ 


(59) 


(60) 


Xk = -k< 

2 ' 


164^X2 


324^X2 


1692 x 2 


329*2x2 


E^El {ojI + <p) + Et^El^ {iol + 9*2) {uj^ + q^) + E^E^, {oj^ + q^) 


+ - 


329-2 (6X'+Xfc+ 60^ I 32g2(6X'+Xfc + y6^ 


El^Et (cc2 + 9-2) 

(iriufe + k) 




169 - 2 x 2 


E4 El (cc2 + 9*2) Et^ (c.2 + 9*2) (aXfe + Yk) 


El 

ao 


+ 


169 - 2 x 2 


169 - 2 x 2 


ElEl {col + 9 - 2 ) ^ ElEl {col + 9 - 2 ) 


169-2(62C^ + ^fc + l4) 


7 \ 2 


869-"X3 


329-2y,'2 


3292yf 


t?^(^^ + 92)(^^Xfc + r6 E4n(cc2+9-*2)(aXfc + n) X4X2(w2+g-2) E^E^ {col + q^) 


+ 


IGq^Yl'^ 


q. W/ X, .W iYC 


E4F,(w2+g2) 
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\iiXkYk-l 


l) 


XkYk+Y^ 

^ i f Miq^Xi _ 4(2Xfc + 3y,0 

2 ]E^El{ujl + q^)^ E^^EtAcol + q^) 


S^iq^Xi 


E-^ 

ao 


E^^Eliul + q^) ElEt{u:l + q^) EtYu{u,l + q^) 

4[n(2eii;" + n')-4ein'' 


EtYk 


Et{u:l + q^){^^Xk+Yk) 


32^iq^Y;^ , 32eig"yf 16a ^ (’ 


XkYk _OTA. 

'a^fc+n 


0 




( 61 ) 


(62) 


with; 


El = k^+u;l + 2Vl + A(,Vl' 

E'^ = k'^ +uj'^ + 214 ' j 


El=k^+u;l + 2Vl + 2^,Wk , 
El = k 
E 


< + 244 , 
ai = + 244 + 2a^fc ■ 


(63) 

(64) 

(65) 

( 66 ) 

(67) 


The meson masses are given by: 


m. = ^214' + 4ai4 


= \/2Vl , 




m, 


a, = sj2V’ + 2a4bfe , 

TOp = \/^ , 


m, 


ai 


= ^244 T2^. 


( 68 ) 
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(70) 

(71) 

(72) 


A.2 Flow equations with U{1)a anomaly 
In the case with U{1)a anomaly, the quantities fk,i in Eq. (1551) read 
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In this case, the meson masses read: 
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(77) 

i?2=fc2+a;2+2i4 + 4fiV," , 

(78) 

El = e+u;l + 2V^ , 

(79) 

Ei^ = e + u;i + 2m, 

(80) 

El = e+u:l + 21^^ + 2fiXfe , 

(81) 

E1 = k‘^ +tol + 2Zk , 

(82) 

El^ = fc2 + u;2 + 2Zk + 2fii4 . 

(83) 


(84) 

1'^" 

II 

t: 

(85) 

moo = V2^/c , 

(86) 

to?7 = \J2Wk + 2fiXfe , 

(87) 

TUp = \/2Zk , 

(88) 

mai = \l 2 Zk + 2fii4 . 
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B Data interpolation 

n data points yj at sites Xj are approximated by a cubic spline /, such that the following expression is 
minimized: 

n 

1% - + (1 -p) 

i=i 

The first term is an error measure, whereas the second a roughness measure. The default value for the 
weights Wj as well as for the weight function A is one. The integration has to be performed over the 
smallest interval containing all data sites, p is a smoothing parameter. This method is used via the 
MATLAB csaps function. 

For both fits in Fig. [TJ all weights are equal to one and p is chosen to be 1 x 10“^. In Fig. [2l we have 
p = 1 X 10“^ and ici = 1 x 10"^ at Xi = 0. 
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